New mechanism for the enhancement of sd dominance in interacting boson models 
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We introduce an exactly solvable model for interacting bosons that extend up to high spin and 
interact through a repulsive pairing force. The model exhibits a phase transition to a state with 
almost complete sd dominance. The repulsive pairing interaction that underlies the model has a 
natural microscopic origin in the Pauli exclusion principle between contituent nucleons. As such, 
repulsive pairing between bosons seems to provide a new mechanism for the enhancement of sd 
dominance, giving further support for the validity of the sd Interacting Boson Model. 
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The Interacting Boson Model (IBM) fl| has been 
highly successful in describing and correlating the col- 
lective properties of medium mass and heavy nuclei 
throughout the periodic table. The model, though phe- 
nomenological in application, is deeply linked to the un- 
derlying nuclear shell model. The s and d bosons of the 
model represent the lowest pair degrees of freedom of 
identical nucleons. The model has a {7(6) group struc- 
ture and three posssible dynamical symmetry limits: 



17(5) D 0(5) D 0(3) 
17(6) D { £17(3) D 0(3) 

0(6) D 0(5) D 0(3) 



(1) 



In each of the three limits, all of which have been re- 
alized to good approximation in nuclei, the hamiltonian 
can be expressed as a linear combination of the Casimir 
operators from the associated group chain. Most nuclei 
do not live at the symmetry limits, however. For a gen- 
eral nucleus, the hamiltonian can be expressed as a linear 
combination of the Casimir operators of all three group 
chains and the nucleus can be represented as a point in- 
side a triangle with the three symmetry limits at the 
vertices M. 

The U(5) chain was initially proposed to describe phe- 
nomenologically vibrational nuclei. Soon after it was pro- 
vided a microscopic interpretation using a mapping pro- 
cedure based on the Generalized Seniority approximation 
M. The other two chains, describing axially deformed 
(£{7(3)) and gamma-soft (0(6)) nuclei, are less well un- 
derstood from a microscopic point of view. 

As noted earlier, the s and d bosons represent the low- 
est pair degrees of freedom of identical nucleons. There 
are of course other pairs as well, but they typically lie 
higher in energy. The key assumption of the IBM is that 
these other pair degrees of freedom can be ignored, ex- 
cept for their renormalized effects on operators in the sd 
subspace. On the other hand, there has been no con- 
vincing microscopic demonstration of this effective de- 
coupling except in vibrational nuclei. 

In this letter, we point out a new mechanism for pro- 
ducing sd dominance in boson models of nuclei. The 
mechanism has at its heart the interaction between alike 



bosons. We will show, by considering a new class of 
exactly-solvable boson models, that the effects of higher- 
spin bosons are suppressed because of the repulsion at 
short range of the boson-boson interaction. Furthermore, 
we claim that a short-range repulsive boson-boson inter- 
action is ubiquitous as it is a direct reflection of the Pauli 
principle at the underlying nucleon level. 

The models we consider involve even angular momen- 
tum bosons extending up to fairly high spin and interact- 
ing via a repulsive boson pairing interaction. As we will 
demonstrate, they can be solved exactly using a method 
introduced by Richardson in the 1960s jf| and recently 
generalized || to confined boson systems. The models 
are closely related to an exactly solvable model proposed 
to describe the transition between the 0(6) and Z7 (5) 
symmetry limits of the IBM [[| . The primary difference 
is that the model of Ref. 6 was restricted to s and d 
bosons, whereas ours accomodate higher-spin bosons as 
well. 

The models we study are based on an £{7(1, 1) group 
algebra, for which the generators are 
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where the operator eJ (et/ m ) creates (destroys) a boson 
with angular momentum Im, and Qi = 21 + 1. We have 
also defined operators ni and A\ which can be used in- 
stead of the generators K. 

In terms of the £{7(1, 1) generators (||), the most gen- 
eral hermitian one- and two-body operator that preserves 
the number of bosons (commutes with the number oper- 
ator) can be written as 



R l =Kf + 2 9 J2{^f 



(3) 
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where the matrices X and Y are at this point arbitrary. 

We will consider the dynamics of an N boson system 
in a Hilbert space that is classified by the product group 
space 

SU(l, l) s x SU(1, l) d x • • • SU(l, l) L , (4) 

where L is the maximum (even) angular momentum per- 
mitted. In this space a model is fully integrable if the 
set oi P — L/2 + 1 global operators R commute with one 
another, viz. 



[Ri,Ri> 







(5) 



We will concentrate in this work on a specific class of 
integrable models. These models correspond to the solu- 
tion of (H) for which Xw = Yw = — j 1 ^ ; , where the rj's 



Vl-Vi' 



are a set of P arbitrary real numbers. The R operators 
of the models can be expressed more compactly as 



Rt (g,r,) = K? 



(6) 



It is straightforward to show that the specific linear 
combination Hpm — 2J2i £ iRi (9: £ ) 01 the R operators 
(^) is precisely the standard pairing model (PM) hamil- 
tonian 



Hpm — £ £lUl 
i 



f£4^+< 
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where c is a constant. 

We now return to the set of R operators (Q) and search 
for their mutual eigenvectors. To do this, we follow 
the procedure introduced by Richardson for the pairing 
model hamiltonian (^). Namely, we consider the ansatz 
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are as many solutions as states in the Hilbert space. 
Thus, the equations (||) define a complete set of eigen- 
states of the form (j|). 

The pair energies e M are always real for boson systems. 
The different states are classified by their configuration 
in the (weak-coupling) limit in which the pair energies e M 
tend to their unperturbed values 2m- In particular, the 
ground state (GS) corresponds to all pair energies in the 
interval 2r] a < e M < 2mi- Excited states can be obtained 
either by breaking pairs or by promoting pairs to higher- 
energy intervals. As an example, the first excited state 
with the same seniority (i.e., the same v) as the GS will 
have M — 1 pair energies in the interval 2r\ s < e M < 2% 
and one pair energy in the interval 2rjd < e M < 2rj g . Com- 
ing back to the set of equations @ , we now see that it 
is a general condition for the common eigenstates of the 
mutually-commuting R operators. Thus, it can be used 
for finding the eigenstates of any linear combination of 
the R operators, not just the linear combination that cor- 
responds to the pure pairing hamiltonian. Once we have 
solved (||) for the pair energies, we can readily obtain the 
eigenvalues of the corresponding R operators from 
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A general pairing hamiltonian can be expressed as a 
linear combination of the R operators as 



H = 2^2 £ lRl(9,v) ■ 



(10) 



Replacing the Ki operators in the Ri operators appearing 
in ( [To|) by the pair operators Ai and the number opera- 
tors ni, and adding and substracting the 5(1, 1); Casimir 
operator, we obtain 



where M = (N — v)/2 is the number of boson 
pairs coupled to zero angular momentum, and \v) = 
\u s , Vd, ■ ■ ■ , vi) is a state with vi unpaired bosons with 
angular momentum I, defined by 

K^\ V )^h\ V ) , K-\v)=0, 

with k\ — v\j2 + fi//4 and v = v\. 

Acting with the R operators on the state (JsJ) and con- 
sidering the conditions that must be satisfied to fulfill 
the eigenvalue equation Ri \*f?) — r\ |^/}, we obtain a set 
of coupled equations for the unknown pair energies e a , 
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This set of equations was first derived by Richardson 
[Ej when considering the eigenstates of the standard pair- 
ing model hamiltonian (m = si). He showed that there 



H = A + ^2 + £ Vw 



A\A V - nfii' 
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where 



a = W ( ^ + 1 E Viin ^ - 1 E Vu ' niQi ' 



(1 



= ei + 2V r u -E V fo' n '' > Vl1 ' =f 



.9 £1 - £v 



2 m - w 



As noted earlier, for E\ = we recover the standard 
PM hamiltonian (|). Restricting @ to I = and 2 
bosons only, we obtain a hamiltonian very similar to that 
of Ref. 6, the only difference being that our hamiltonain 
has an additional monopole-monopole interaction term 
that can be traced to the term involving Yu> in (||). 
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FIG. 1. Occupation probabilities for the ground state of 
a system of M — 5 boson pairs and maximum angular mo- 
mentum L — 12 as a function of the interaction strength g. 
The upper graph shows the sum of occupation probabilities 
(depletion) for high-spin bosons (I > 2) while the lower graph 
gives the occupation probabilities for s and d bosons. 

An important feature of our exactly solvable and inte- 
grable hamiltonian ([ll]) is that it is not restricted to an 
sd subspace but also permits higher angular momentum 
bosons. We will now make use of this property to study 
the effect of the inclusion of high-spin bosons on the sd 
subspace dynamics. 

Since there is virtually no information available on the 
coupling of the usual IBM space to higher-spin bosons, 
we will consider two sample scenarios. In the first (de- 
noted I), we assume a constant pairing interaction, as 
arises for r\i = We also assume linear single-boson 
energies (e; = I). 

As is well known, a pure pairing interaction produces 
enhanced pair scattering to high-spin states. This is pre- 
sumably a consequence of the extreme short-range nature 
of a pairing force, and is reflected in the fact that the ef- 
fective pairing matrix element that connects a level with 
angular momentum I to one with angular momentum I' is 
Vi e J f = Viv % /(2Z + l)(2/' + l) . This anomalous feature 
of a boson pairing force bears some similarity to the well- 
known pathology of a delta interaction in fermion theo- 
ries. There, the unphysical behavior is typically overcome 
by introducing a cutoff in energy. In our boson context, 
we can likewise modify the interaction to get rid of its 
undesirable high-^ behavior. In particular, if we select 
the rj's according to r]i — I 2 , the effective pairing matrix 
elements would scale as 



V - 



eff 9y /(2l + l)(2l' + l) 



2(1 + 1') 



and not produce enhanced scattering to high-spin pairs. 
In the second scenario (denoted II), we will make this 



assumption for the r/'s and again assume linear single- 
boson energies (e/ = I). While we feel that model II 
is more realistic than model I, we will for completeness 
present results for both. 

For both models I and II, as well as for any other gen- 
eral pairing model, the occupation probabilities can be 
readily obtained from the R operators according to 

Using the Hellmann-Feynman theorem, we can replace 
the R operators by their eigenvalues to obtain an equa- 
tion for the occupation numbers in terms of the deriva- 
tives of the pair energies with respect to the pairing 
strength g. Furthermore, taking the derivative of the 
Richardson equations (^) with respect to g leads to a set 
of linear equations for the pair energy derivatives. In Fig. 
1, we show the occupation probabilities for the GS of a 
system of 10 bosons and an angular momentum cutoff of 
L = 12 as a function of g. In the lower part of the figure 
we show the occupation probabilities for s and d bosons, 
while in the upper part we show the summed occupation 
probabilities of those bosons with I > 2, which we call 
the depletion of the IBM space. In both parts of the fig- 
ure, the dashed line refers to model I and the solid line 
to model II. 

It is apparent from the figure that the class of models 
we have considered displays a quantum phase transition 
from an s boson condensate (spherical) to a mixed state 
of s and d bosons ("deformed"). The transition is soft- 
ened by the finite number of bosons so that there is still 
some occupation probability of high-spin bosons. Two 
important features should be noted in these results. One 
is that the depletion increases in the spherical phase up 
to the phase transition and then stays almost constant 
in the deformed phase. The other is related to the choice 
of the model interaction. In the first scenario with the 
PM hamiltonian and linear single-boson energies the de- 
pletion is small but non-negligible; in the more physical 
second scenario the depletion is negligible and sd domi- 
nance is almost complete. 

The fact that these models all lead to a very small non- 
sd content in their ground states is the key conclusion of 
this paper. It suggests that repulsive pairing is indeed a 
robust mechanism for reinforcing sd dominance in boson 
models of nuclei. 

The above analysis was for a system of 10 bosons. An 
important issue not addressed by those results concerns 
the dependence on the number of bosons. How many 
bosons are needed for the repulsive pairing interaction to 
be effective in suppressing the high-spin content of the 
low-lying states? 

This question is addressed in Fig. 2, where we show 
the s and d occupations and the IBM depletion for the 
GS as a function of the number of boson pairs M for the 
specific value of g = 0.5. As is evident from Fig. 1, this 
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choice of the pairing strength places the system within 
the deformed phase, where there is mixing of different 
bosons in the ground state. 
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FIG. 2. GS occupation probabilities for a system of bosons 
and maximum angular momentum L = 12 and pairing 
strength g = 0.5 as a function of the number of boson pairs 
M. The upper graph shows the sum of occupation probabil- 
ities for high-spin bosons (I > 2) while the lower graph gives 
the occupation probabilities for s and d bosons. 

The first point to note is that for small enough bo- 
son number, the depletion from the sd subspace actually 
increases with increasing boson number. What is impor- 
tant, however, is that for both hamiltonians a critical 
boson number is soon reached after which the depletion 
decreases and sd dominance gradually becomes better 
and better. The precise number of bosons required for 
this critical behavior to set in depends on the hamilto- 
nian. It is 12 for the standard pairing interaction and 6 
for the more realistic modified pairing interaction. 

It is interesting to ask why in the presence of a repul- 
sive pairing interaction there is a phase transition to a 
state of essentially pure sd content. 

As mentioned before, the GS solution of (||) for repul- 
sive pairing corresponds to all pair energies in the interval 
2rj s < e M < 2rjd. This means that the pair operators in 
the Richardson ansatz (|J) have one phase for s boson 
pairs and the opposite phase for all other boson pairs. It 
is only when this set of phase relations is satisfied that 
the system can gain in energy from a repulsive pairing 
interaction. On the other hand, for two given boson de- 
grees of freedom, the only way they can take advantage 
of a repulsive pairing force is by having opposite phases. 
Putting these two facts together, we see that only one 
of the non-s boson degrees of freedom can correlate with 
the s boson and produce a gain in energy in the presence 
of a repulsive pairing interaction. Clearly, this will be 
the d boson since it has the lowest energy. The picture 



that emerges is that the d boson mixes with the s boson 
once the pairing interaction is strong enough to overcome 
the difference in single-boson energies, thus explaining 
the phase transition from SU(b) to 0(6) in the IBM as 
a function of the pairing strength. Higher-spin bosons 
do not lead to significant further correlations, since they 
cannot have opposite phases to both the s and d bosons 
that make up the 0(6) wave function. 

At this point it is worth expanding briefly on the com- 
ment made earlier that a short-range repulsive interac- 
tion between composite objects is a natural consequence 
of the Pauli principle at the constituent level. That 
this is the case is well known both in molecular physics, 
where the short-range repulsion between the two hydro- 
gen atoms in an H2 molecule is the result of electron ex- 
change, and in nuclear physics, where the nuclear force 
has a strong short-range repulsion due to quark exchange. 
It also arises naturally when boson mapping methods are 
applied to systems with spatial two-body correlations . 

In this work, we have discussed a class of boson mod- 
els involving a generalized repulsive pairing interaction 
for bosons with arbitrary even angular momenta. We 
have shown that these models are exactly solvable using 
a method originally developed by Richardson for more 
traditional pairing hamiltonians. We have applied this 
method to a system in which the boson degrees of free- 
dom extend to a cutoff angular momentum significantly 
larger than 2. Nevertheless, we find that the ground state 
of the system is almost completely dominated by s and 
d bosons, the effect of higher-spin bosons being strongly 
suppressed by the repulsive pairing interaction. Consid- 
ering that a repulsive pairing interaction between bosons 
is a natural consequence of the Pauli principle at the 
nucleon level, these results provide further microscopic 
support for the validity of the interacting boson model 
in nuclear structure. 
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